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A macroscopic equation of motion for the magnetization of a ferromagnet at elevated temperatures 
should contain both transverse and longitudinal relaxation terms and interpolate between Landau- 
Lifshitz equation at low temperatures and the Bloch equation at high temperatures. It is shown 
that for the classical model where spin-bath interactions are described by stochastic Langevin fields 
and spin-spin interactions are treated within the mean-field approximation (MFA), such a "Landau- 
Lifshitz-Bloch" (LLB) equation can be derived exactly from the Fokker-Planck equation, if the ex- 
ternal conditions change slowly enough. For weakly anisotropic ferromagnets within the MFA the 
LLB equation can be written in a macroscopic form based on the free-energy functional interpolat- 
ing between the Landau free energy near Tc and the "micromagnetic" free energy, which neglects 
changes of the magnetization magnitude |M|, at low temperatures. [80163-1829(97)03905-2] 
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I. INTRODUCTION 

The famous Landau-Lifshitz equation,0 which is the 
basis of innumerable investigations of magnetically or- 
dered materials, considers magnetization as a vector of 
fixed length and ignores its longitudinal relaxation. Such 
an approach is obviously unsatisfactory at elevated tem- 
peratures since magnetization is an average over some 
distribution function and its magnitude can change. Al- 
ternatively, semiphenomenological "soft-spin" equations 
of motion for the spin density allowing for the longitudi- 
nal relaxation and for the influence of the bath described 
by stochastic Langevin terms ai£_kjiown in the theory 
of dynamic critical phenomena.0'&a A phenomenologi- 
cal deterministic equation of motion for the magneti- 
zation of magnetically ordered materials with the lon- 
gitudinal relaxation terms, which is a direct generaliza- 
tion of the Landau-Lifshitz equation, was formulated by 
Bar'yakhtard'H and applied to the domain- wall dynam- 
ics at elevated temperatures.!] The Bar'yakhtar equation 
was conceived for the temperature range below the Curie 
point Tc; the theory does not answer what happens 
with phenomenological relaxation terms above Tc and 
whether the Bloch equation used in the theory of EPR 
and NMR can be recovered in this region. 

The simplest nontrivial model, for which the prob- 
lem of finding an equation of motion for magnetiza- 
tion in the whole range of temperatures can be formu- 
lated, is a semiphenomenological model considering an 
isolated classical spin interacting with the bath modeled 
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by stochastic Langevin fields. The spin-spin interactions 
in this model, which lead to the ferromagnetism, can be 
taken into account on the next stage on the mean-field 
level. Dynamics of such a spin is described by the Fokker- 
Planck equation (FPE) , which can be solved analytically 
only in limiting cases, in particular, of low and high tem- 
peratures. Reduction o£-|the FPE using the modeling of 
the distribution functions (the accuracy of this procedure 
was shown to be about 7% in most situations) has led to 
the closed equation of motion for magnetization interpo- 
lating between the Landau-Lifshitz and Bloch equations 
at low and high temperatures — the so-called "Landau- 
Lifshitz-Bloch" (LLB) equation. The LLB equation was 
also derived for a quantum spin system interacting with 
a bathtl by the reduction of the density-matrix equation 
with the method similar to that used in the classical case. 
A kind of LLB equation taking into, account the spin- 
spin relaxation was obtained by PlefkacjEil for a quantum 
model with long-range "spin-block" interactions. 

The coefficients in the relaxation terms of such a gen- 
eral LLB equation are nonlinear functions of magneti- 
zation itself; the only application of this equation up to 
now is that to the calculation of the nonlinear mobility pi 
domain walls (DW) in rare-earth (RE) fcrrites garnets,Ej 
where the strongly thermally disordered spins of the RE 
sublattice do not interact with each other and are sub- 
ject to only the combined influence of the external field 
and the molecular field acting from the iron sublattice. 
For the simplest one-sublattice weakly anisotropic fer- 
romagnetic model below Tc the dominant term in the 
molecular field is the homogeneous exchange, so that the 
directions of the molecular field and magnetization nearly 
coincide. In this case the general LLB equation sim- 
plifies to its particular form similar to the Bar'yakhtar 
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equation. The latter was applied in Refs. p| to 
calculate the domain-wall mobility in uniaxial ferromag- 
nets in the whole temperature range and, in particular, 
near the phase transition from Bloch to linear (Ising- 
like) walls at snmc Tb < Tc predicted by Bulaevskii 
and Ginzburg.Ej As this second-order phase transition 
is accompanied by changing the roles of transverse and 
longitudinal relaxation processes in the DW|-dMaaniics, 
the DW mobility has a deep minimum atTsMBElThis 
minimum, and thus the DW phase transition, was re- 
cently observed in dynamie-ipjficeptibility experiments on 
Ba and Sr hexaferrites.llaolla 

An important dynamical scenario is that when the rate 
of changing of magnetization (or of its spatial distribu- 
tion), which can be controlled by an external influence, 
is slow in comparison to the spin-relaxation rate. This 
small parameter makes it possible to solve the Fokker- 
Planck equation exactly without making assumptions 
about the form of the distribution function. For ex- 
ample, calculation of the low-frequency imaginary part 
of the longitudinal susceptibility leads to the exact an- 
alytical expression for the integral relaxation time Tint, 
which is defined as the area under the magnetization 
relaxation curve afi|ejpa]H-ab|rupt infinitesimal change of 
the magnetic field.&y'BEl'E3 The quantity Tint describes, 
in particular, the thermoactivation escape rate of fine 
ferromagnetic particles with a uniaxial anisotropy over 
the potential barrier, which is-,valid, in contrast to the 
well-known Brown's solution,Bj in the whole tempera- 
ture range. Such a situation is also characteristic for the 
dynamics of domain walls, whose velocity depends on 
the amplitude of the driving field and can be kept what- 
ever small. In this case the Fokker-Planck equation can 
be solved exactly, which leads to the exact form of the 
LLB equation, if the spin-spin interactions are considered 
within the mean-field approximation (MFA). Derivation 
of this exact "slow" form of the LLB equation is the main 
purpose of this article. 

The main part of the paper is organized as follows. In 
Sec. H the Fokker-Planck equation for a classical spin, 
its low- and high-temperature solutions, and the approx- 
imate reduction of the FPE to the Landau-Lifshitz-Bloch 



equation is outlined. In Sec. Ill the FPE is exactly solved 
in the slow-motion case and the slow LLB equation is 
derived. In Sec. IV the simplified form of the latter for 
ferromagnets below and near Tc is worked out. In Sec. 
further possible applications of the method and some 
unsolved problems are discussed. 



II. THE FOKKER-PLANCK AND LLB 
EQUATIONS 



In the case of a weak coupling with the bath the dynamics 
of the vector s can be described with the help of the 
stochastic Landau-Lifshitz equation 



s = 7[s X (H + C)] - 7A[s X [s X H]; 



(2.1) 



with A <C 1, where correlators of the a, (3 = x,y, z com- 
ponents of the Langevin field Cit) ^re given by 
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(2.2) 



The Fokker-Planck equation corresponding to Eq. 
( p.lD is formulated for the distribution function /(N, t) = 
((5(N — s{t))) on the sphere |N| — 1, where the average 
is taken over the realizati ons of C- Differentiating / over 
t with the use of Eq. (2.1) and calculating the right 
part with the methods of stochastic theory (see the Ap- 
pendix) , one comes to the Fokker-Planck equatiorfl 
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One can easily see that the distribution function 
/o(N) cx exp[-7i(N)/r], His) = -^qHs 



(2.3) 



(2.4) 



satisfies (2.3) at an equilibrium. For the small coupling 
to the bath, A ^ 1, Eq. (2.3) coincides with the Fokker- 
Planck equation derived by Brown.cj 

The equation of motion for the spin polarization (the 
first moment of the distribution function) 



m = (s) = /d^AfN/(N,t) 



(2.5) 



of an ass embly of magnetic atoms can be derived from 
Eq. (2.3) and has the form 



m = 7[m X H] - AArm - 7A([s X [s X H]]) (2.6) 



[cf. Eq. (2.1)], where Ajv is the characteristic diffusional 
relaxation rate or, for the thermoactivation escape prob- 
lem, the Neel attempt frequency given by 



An = 



-1 _ 



N 



= 27Ar/A*o. 



(2.7) 



It can be seen that Eq. (2.6) is not closed but coupled to 
the second moments of the distributio n fu nction, (siSj), 
in its last term. The behavior of Eq. (2.6) is determined 
by the reduced field given by 



^0 = MoH/r. 



(2.8) 



We shall describe a magnetic atom as a classical spin 
vector s of a unit length. The magnetic and mechanical 
moments of the atom are given by /x = ^qS and L = 
fiQs/j, where 7 — ge/{2mec) is the gyromagnetic ratio. 



For ^0 3> 1 (low temperatures) the second term in Eq. 
( p.6D can be neglected and the last term decouples for 
distribution functions localized about some direction: 
(siSj) = mirrij. In this case the Landau-Lifshitz equa- 
tion of the type (2.1) for m without the stochastic field 
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^ is recovered. In the high -temperature case, ^ the 
second term of Eq. (2_^) dominates over the last one, 
which can be neglected. Here one gets the equation of 
motion for m with the Bloch relaxation term. 

In the intermedi ate region, ~ 1; where the first- 
moment equation (2^) is not close d, t he resonance and 
relaxational behavior of the FPE (2.3) is not described 
by Lorentz and Debye curves., and the deviations from 
the latter reach 7% at 3.1a Neglecting these features, 
one can obtain an isolated equation of motion for the spin 
polarization of an assembly of magnetic atoms choosing 
the distribution function in a formB 



/(N,t) 



exp[g(t)N] 



Z ^Att 



sinh^ 



(2.9) 



cf. Eqs. (|2^ and ( p| ) 



where ^(t) is chosen so that the 
first moment equation ( 2.6| ) is satisfied. Calculating the. 
terms of Eq. (2.6) with the help of Eq. (2.9), one arrivesO 
at the LLB equation for the nonequilibrium reduced field 

m 



^ = 7K X H] - Ti ( 1 - ^ ) ^ 



(2.10) 



with the longitudinal and transverse relaxation rates 
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(2.11) 



where Aat is given by Eq. (p7|), B{^) = cothf ~ ^/C 
is the Langevin function and B'{^) = dB{£^)/dS^. The 
asymptotic forms of Fi and r2 are given by 




and 



A. . . , 



(2.13) 



The relaxation rates of such a type also appear as a result 
of calculation of the high-frequency longitudinal suscep- 
tibility and the far-from-resonance transverse one.u The 
quantity Ti is also proportional to the "effective eigen- 
value" Aef of Ref_0. One can see that the equilibrium 
solution of Eq. (2.1C) is ^ = ^o- The nonequilibrium spin 



polarization m is given by 
m = mi/i, 



(2.14) 



The LLB equation for |, Eq. ( |2.10| ), can be written in 
the alternative equivalent form 



i ^ 7[l X H] - r,(| - ^o) - (r^ - ro ^^^l^ 

(2.15) 

Here it can be seen that in the high-temperature region, 
C^o < 1, where B{£,) ^ ^/3 and Li ^ Ls = Aat, the 
Landau-Lifshitz double-vector-product relaxation term 
becomes small and the Bloch equation is recovered. On 
the other hand, at low temperatures, when ^,^o ^ li 
the magnitude of the vector m in Eq. ( 2.14 ) saturates in 
most situations at m — B(£) = 1, and the longitudinal 
relaxation term in Eq. (2.10) no longer plays a role. Here 
the usual Landau-Lifshitz equation is recovered. Using 
Eq. (2.14) one can derive the LLB equation for the spin 
polarization m itself. The result can be written as 



rh = 7[m x H] — Aat 1 



7A 1 - 



m X [m X H]' 



(2.16) 



[cf. Eqs. {2A) and ( |2.6| )], where ^ = is determined 
implicitly by the relation m — B{^). Note that here 
at low temperatures, ^ ^ 1, the coefficient before the 
transverse relaxation term goes to 7A, whereas the lon- 
gitudinal one is nonessential, if m is saturated. At high 
temperatures, ^ <IC 1, the relaxation term in Eq. (2. If) 
acquires the Bloch form A7v(ni — mo) with mo = $o/3 
[see also Eq. (2.15)]. The quantum generalization of the 
classical LLB equation written above was given in Ref. 
^. The latter was applied in Ref. |l^ to study the nonlin- 
ear dynamics of the RE sublattice of rare-earth ferrites 
garnets near the magnetization compensation point. 

For small deviations from equilibrium, where ^ ^ 
and, accordingly, m ^ mp = -B(Co)Co/'?o, one can put 
the LLB e quation (2.16) [or, more conveniently, directly 



(2.12) (I2.IOD ] into the form 



m = 7[m X H] — Fi (^1 - 
[m X [m X mo 



mmo\ 



m 



F2 



(2.17) 



where the relaxation frequencies Fi and F2 are functions 
of ^0- A kind of LLB equation similar to Eq. (2.17) 
was obtained by Gekht et al.f3 who assumed, for the 
calculation of the lin ear transverse dynamic susceptibil- 
ity, instead of Eq. ( ^.9| ) a distribution function of the 
form /(N,0 = /o(N)[l + a{t)N], where fo is given by 
Eq. (2^) and ct corresponds to ^ — ^0 in our notations. 
Although Gekht et al. claimed that "the single-moment 
approximation is permissible for small deviations from 
equilibrium," Eq. (2.17) is in fact only approximate, as 
well as the more general Eq. ( 2.16| ). The latter, in con- 
trast, can be applied and has a rather good accuracy in 
situations where deviations from equilibrium are large, 
as was checked in R ef. 1 1 . In Sec. HI we will consider the 
solution of the FPE (|2.^ for slowly varying field H(<). In 
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this case the deviations from the instantaneous equiHb- 
rium state are small and the FPE can be solved exactly 
without assumptions about the form of the distribution 
hmction f{N,t). 



III. THE "SLOW" LLB EQUATION 

If the magnetic field H slowly changes its magnitude 
and direction, the solution of the Fokker-Planck equation 
(2.3) slightly deviates from the instantaneous equilibrium 
one and can be searched for in the form 



Q<1, (3.1) 



where ^o(^) = MoH(i)/T. The correction function 
(5(N,t) oc |H| and, additionally, it depends slowly on 
time, so that Q cx | Hp . Neglecting this small term, one 
obtains from Eq. (2.3) the equation for Q having the 
form 
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ro(mo - N)^o, mo = B{^q) 



H 



(3.2) 



where tq = ^^/{■^T). One can see that in leading order 
the correction (3(N, t) is determined by the instantaneous 
values of the magnetic field H(i) and its first derivative 
H. The right-hand part of this equation can be separated 
into the terms describing the temporal changes of the 
magnitude and of the direction of H as 



(mo - N)^o = N[^o X f2] + mo 
where 

= [lo X iol/Co 



N|o 

Co 



io (3.3) 



(3.4) 



is the precession frequency of the vector ^o- In th e sp her- 
ical coordinate system with z axis along ^o Eq- ( ) for 
(3(x, (p), where x = cos0, takes on the form 
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Q 



— To^o \/l — [Vly COS ip — 

+ Ta{x ~ mo)Co, 



flxSmcp] (3-5) 



where Vl^ and Qy are x and y components of the vector 
ft. 

The solution of the linear differential equation ( |3.5| ) is 
a sum of two contributions induced by the transverse and 
longitudinal inhomogeneous terms: Q = Qj_ + Q\\. Using 
the substitution 



Q± = Qx cos (p 



(3.6) 



one comes to the equation 
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dx 



dx 



0^ 



(3.7) 



where fi+ = fli + iD,2- This equation cannot in general 
be solved analytically, but the latter is possible in the 
typical case of the weak coupling to the bath, A ^ 1. 
For A/^o 1 one can easily find the solution iteratively, 
which yields 



Q+ ~ Tafl+\/T~~x^ 



1 



iX 
1^ 



(2 + ^ox) 



(3.8) 



On the other hand, in the high-temperature region, where 
f -C 1, one can neglec t £o in the round brackets in Eq. 
( |3.7|) , after which Eq. (3/7) can be analyticaly solved to 
yield 



To 



: 1 + 2zA/eo 
'l + (2A/eo)2 



(3.9) 



These two solutions overlap in the region A ^ ^o *C 1, 
and thus they can be sewn together in the whole range 
of temperatures into the formula, which can be obtained 
by replacing the numerator of the fraction in (|3.9| ) by 
l + (zA/eo) (2 + ^02^). 

The equation for Q\\{x) can be written as 



;^ + eo)(l-x2)^ 
dx / dx 



= A]^Hx - m„)io- (3.10) 



It can be solved in two steps with the help of the sub- 
stitution P{x) = (1 — a;2)d(5||/dx. First, integrating Eq. 
( |3.10 ) one gets 
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Then, is given by 

Qii(^) = 



coth ^0 



sinh ^0 



dx' 



1 



„/2 



P{x) + C, 



(3.11) 



(3.12) 



_1 -L — X'- 

where the constant C is determined from the normaliza- 



tion condition for the distribution function (3.1). 

Now, the function (5(N,t) having been determi ned , 
one c an c alculate the spin polarization m using Eqs. ( |2.5| ) 
and (3.1). Returning to vector designations, one comes 
to the result 



5(^0 




i^B' HH ^ H 
H 



[H X H] Fa [H X [H X H]; 



(3.13) 
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where r2 is the transverse relaxation rate given by Eq. 
( ^.11 ) and Fi^int is the inverse of the integral longitudinal 
relaxation time Tint, 



1 



Tint 



int 



AAr^osinhCo-B'(Co) 



dx 



coth^o 



sinh ^0 



1 -a;2 

(3.14) 



which is determined as the area under the magnetization 
relaxation curve after an abrup|t|iHfinitesimal chang e of 
the longitudinal magnetic field.Eli23 Equation (3.13) de- 
scribes the lagging of the spin polarization m from its 
quasiequilibrium value mo(t) of Eq. ( |3.2| ), which is de- 
termined by the small deri vative H(i). The asymptotic 
forms of Fi^iiit in Eq. ( 3.14| ) read 
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Tl,: 
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Aa'Co 1 



Co 



Co » 1- 



(3.15) 



Comparing Eqs. ( 3.15| ) and (2.12) one can see that Fi > 
Fi,int- The relative deviation d = Fi/Fi — 1 attains a 
vaiue S « 0.07 at « 3B 

The next problem is to write dow n the equation of mo- 
tion for m, which has the solution ( 3.15 ). It is especially 
important if the spin-spin interactions are taken into ac- 
count within the MFA (sec the next section). In this case 
H is replaced by the m olecular field H^^-^ containing m 



itself, and Eq. ( |3.13| ) is in fact a differential equation 
for rh, which shoul d be still simplified. It can be done 
differentiating Eq. ( |3.13 ) over time and neglecting terms 
of or der coming from the correction terms with H in 
Eq. ( |3.13D . This leads to 



,(HH)H .^JHxfHxHll , 



H3 



Now H in this relation shou ld be expressed through m 
with the help of Eq. (3.13), which after some vector 
algebra leads to the "slow" LLB equation 



m = 7[m X H] - Fi.int [1 
[m X [m X mo]] 



mmo 



(3.17) 



where mo is g iven by Eq. (3.2) and which is the refine- 
ment of Eq. ( 2.17 ) in the s low-motion situation. The 
quantities Fi of Eq. ( ^.ll[ ) and Fijnt of Eq. ( |3.14| ) 
have the same leading high- and low-temperature asymp- 
totes, and, as was said above, they differ by no more 
than 7% in the whole range of temperatures. The same 
order of magnitude also characterizes the difference be- 
tween the Debye one-relaxator form of the longitudinal 
dynamic susceptibility x\\{^) following from Eq. ( |2.17 ) 
and the actual form of X||(w) following from the solition 



of the exact Fokker-Planck equation (2.3) at intermediate 
temperatures.l3 It sh ould be noted that in the fast-m otion 
situations equation (2.17) is better than Eq. ( |3.17 ), since 
it yields the exact leading (imaginary) term of the high- 
frequency expansion of Xj|(Li;).El 



IV. LLB EQUATION FOR FERROMAGNETS 

For definiteness we consider the classical ferromagnetic 
model with the biaxially anisotropic exchange interaction 



-/io ^ Hi 



rjySyiSyj 



(4.1) 



where yy^. < rjy < 1 are the anisotropy coefficients. The 
dynamics of this model interacting with the bath is de- 
scribed by the stochastic Landau-Lifshitz equation 



s» = 7[si X (Hi,tot + Ci)] - 7A[si X [si x H;, 



(4.2) 



[cf . Eq. ( p.l[ )] , where Ci are postulated to be uncorrelated 
on different lattice sites, and 



H, 



1 dH 

Mo dsi 

A^o 



(4.3) 



is the total field acting on a given spin at the site i, which 
depends on the orientation of spins on the neighboring 



sites j. In Eq. ( |4.3| ) Sq 



X, y, z, and are 



the orts of the Descarte coordinate system. 

The Fokker-Planck equation for the distribution func- 
tion 



AT 



/tot({N,},t) = (n'5(N,:-s,(t)) 



(4.4) 



of the whole system cons istin g of J\f spins can be derived 
in the same way as Eq. ( |2.3| ) and has the form 

dfu 



dt 



h ^ [N. x H,,tot] - 7A[N, X [N, x H.^tot]] 



Mo 



N,: X 



N,: X 



d 
9N" 



fu 



0. 



(4.5) 



One can check that the static solution of this equation is 
/tot,o({Na) cx exp[-H{{N,})/T] (4.6) 



where H is given by Eq. (4.1). Solving Eq. ( |4.5| ) is a 
formidable task that goes beyond the scope of this paper. 
It is in any case not simpler than calculating averages 
with the distribution function (4.6) at an equilibrium and 
requires application of some kind of many-body pertur- 
bation theory, as the diagram technique for classical spin 
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systems (see, e.g., Ref. ^5|), which has proved to be rather 
efficient for description of their static properties. Here we 
resort to the mean field approximation with respect to 
spin-spin interactions, which means, however, dropping 
their contribution into the relaxation rates. In MFA the 



distribution function of the system (^^) is multiplica- 
tive, and one can use the distribution functions fi for 
each spin on the site i, which satisfy the Fokker-Planck 
equatio n dpi ) with H Hf^^^, where Hf^^^ is given 
by Eq. ( |4.3| ) with the replacement S; mi = (s^). Solu- 
tion of such mean-field FPE's similar to that of Sec. Oor 



Sec. Ill leads to the set of coupled LLB equations for nii, 
i = 1, 2, ... ,7V o f the type ( p. 16 ) in a general nonlinear 
situation or Eq. ( 3.17 ) for slow motions. The static solu- 
tion of these LLB equations satisfies the inhomogeneous 
Curie- Weiss equation. 



Mo 



H 



MFA 



T 



(4.7) 



which describes within the MFA both the homogeneous 
state and such configurations as domain walls with ac- 
count of thermal effects (see, e.g., Ref. and references 
therein) . 

For the most of ferromagnetic substances the small- 
anisotropy case, i.e., rj'^ y = 1 — ij^^y ^ 1, is realized. In 
this case the spatial inhomogeneity of magnetization at 
a distance of the lattice spacing is small, and one can 
use the continuous approximation. For H^^^ the latter 
means 



H 



MFA 



He = — m, 

Mo 



Jn 



H^ff = H + ^[aAm-77>,-r;;m^], (4.8) 
Mo 

where Jq is the zero Fourier component of the exchange 
interaction, A is the Laplace operator, and a is a lattice- 
dependent constant (for the simple cubic lattice a — 
Oq/G, where oq is the lattice spacing). The most im- 
portant for ferromagnets is the case of the strong homo- 
geneous exchange field, |H£;| ^ IH^ffl, which is reahzed 
below Tc = ^Joi where there is a spontaneous magne- 
tization, and also in the region just above Tc, where 
the longitudinal susceptibility is large. As in this case 
the external field H(i) that can drive the system off the 
equilibrium is a relatively small quantity, one can use 
Eq. (2.17) [or, for slow motions, Eq. (3.17)] and expand 
mo = B{pfioH^'^^)Ii^^^/H^^^, where /? = 1/T, up to 
the first order in H^g. This leads to the equation 



rh = 7[m x H^^] — 7A 



1-B/m mH'^ 



off 



-7A5 



M0/5S' 
[m X [m X H'^g]] 



m 



where B = B['m[3Jo), 



Ai = 2A 



T 
Jo' 



A 1 



T 
Jo 



(4.9) 



(4.10) 



if Eq. (2.17) was used, and the sam e w ith Ai 
AiFiant/Li for the "slow" LLB equation (|t^). The dif- 
ference 1 — _B/m in Eq. ( |4.9| ) is a small quantity propor- 
tional to the deviation from the equilibrium. It can be 
further simplified to 



1 - B/m 
Mo/35' 



1 

Jo 
Mo 



T <Tc 
|e| « 1, 



(4.11) 



where e = 1 — T /Tc, me is the equilibrium spin polariza- 
tion satisfying rrie — B{me(3Jo), and 



XII 



dm /io B'(3Jo 
Im ^ Jo 1 - B'pjo 



(4.12) 



is the spin polarization susceptibility, calculated for m — 



me- 



Using B{0 = U- 



and 



jc near 



Tc, one can check that the two expressions in Eq. (4.11) 
overlap in this region. 



The last step is to rewrite (4.9) for the macroscopic 
magnetization, M = fiQin/vQ, where vq is the unit-cell 
volume. This leads to the final result 



M = 7[M X Hcff] + Li 



(MHeff)M 



M2 



[M X [M X Heff]] 
' Af2 



(4.13) 



where Li and L2 are the longitudinal and transverse ki- 
netic coefficients. 



(4.14) 



cki and ai are the corresponding Gilbert damping pa- 
rameters, and the effective field Heff is given by 



H 



1 

7, 



AM 



1 



1 fAP 



Xx 

1 1 M 



1 

Xy 



[cf. Eq. (|4^)]. In Eq. ( [4. 151) 



aJo 
Wd' 



Mo 



(4.15) 



(4.16) 



Qd and Wd are the characteristic dipolar wave number 
and dipolar energy, a Jo is the second moment of the 
exchange interaction, and the susceptibilities are given 
by 



XII 



Wd B'pJo 
Jo 1 - B'pjo 



Xx, 



Wd 



1 



Jo 1 ^ 11x,y 



(4.17) 



The effective field U^s of Eq. ( |4T5| ) can be written as a 
variational derivative 



Heff(r) = 



5F 
<5M(r) ' 



(4.18) 
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where F is the MFA free energy of a ferromagnet, 



dr I -HM 



1 



2 (VM)2 



1 

2x. 



-Mi 



— Af2 

2X. ^ 



8M|x|| 



(4.19) 



(VM)2 = (VM^)2 + (VMj^)2 + (VM^)2, and Fq is the 
equiUbrium free energy in the absence of anisotropy and 
magnetic field. The direct derivation of this free en- 
ergy from the mean field theory is tricky and will be 



presented elsewhere. Equation 1 ) provides a link be- 
tween the "micromagnetics" jatA which ignores changes 
of the magnetization magnitude |M|, and the Landau 
theory of phase transitions ^lSEj which is a limiting form 
of the MFA pretending to be valid only in the vicinity of 
Tc where the order parameter M(r) is small. In fact, for 
weakly anisotropic systems in a magnetic field smaller 
than the homogeneous exchange field He, the actual 
small quantity, which remains small in the whole tem- 
perature range, is not Af^(r), but rather the difference 
(r) — Mg , where Mg is the equilibrium magnetization 
in the absence of anisotropy and magnetic field. Since in 
the MFA near Tc one has oc x|^^ cx e = 1 - T/Tc, 

the last term of Eq. (4.19) takes on the Landau form 
AM"^ + BM^ with A = -e^g , and Aq, B = const. This 
shows, further, that Eq. ( 4.19| ) can be continued into the 
region T > Tg as the usual Landau theory. The free en- 
ergy Eq. (4.19) can be brought into the "micromagnetic" 
form by introducing the magnetization direction vector 
ly = M/M. One can then identify 



K 



(4.20) 



where Kx.y are the anisotropy constants. 



V. DISCUSSION 

In this paper several forms of the Landau-Lifshitz- 
Bloch (LLB) equation of motion for a single classical spin 
interacting with the bath as well as for classical ferromag- 
nets within the MFA have been obtained. These LLB 
equations are applicable for all temperatures and contain 
both transverse and longitudinal relaxation terms. The 
nonlinear response of a single spin to the arbitrary chang- 
ing magnetic field H(t) is the most accurately described 
by the nonlinear LLB equation ( 2.16| ). Fo r slo wly vary- 
ing H(i) the exact "slow" LLB equation (3.17) contain- 
ing the integral longitudinal relaxation time Eq. ( 3.14 ) 
can be used. This case is the most important one for 
the domain- wall dynamics. For ferromagnets within the 
MFA the magnetic field H in the LLB equati on s hould 
be replaced by H^^^, which is given by Eq. (4.3) with 
Si \n.i = (si) in a general case or by Eq. (4.8), if the 
continuous approximation is applicable. If, additionally, 



in Eq. (4.8) the homogeneous exchange field H^; is dom- 
inant, which is typical for ferromagnets below and near 
above Tc, t he L LB eq uation reduces to t he fo rm (4.13) 



with Eqs. (tl.lBI) and ( [4.191) . Equation ( [1.13D could be 
also written, without specifying the form of kineti c coef- 
ficients Li 2 and that of the free energy Eq. ( 4.19 ), from 
general arguments. Itj/gry close to the phenomenological 
Bar'yakhtar equation,Efu which contains an additional re- 
laxation term proportional to Aifeff. This term, whose 
microscopic origin is the spin-spin interaction or the cor- 
relation of the Langevin fieldSnCi in Eq. (4.2) on different 
lattice sites i ^ j , was shown,Q however, to yield a contri- 
bution into the domain- wall dynamics, which is negligible 
in comp ariso n to that of the longitudinal relaxation term 
in Eq. ( [4.13|) . 

The qu antum generalization of the nonlinear LLB 
equation ( 2.16| ), which contains additional relaxation 
terms of a different symmetry, was derived in Ref. ^ by 
the approximate solution of the density matrix equation 
for a single spin interacting with an idealized phonon 
bath, which is based on choosing the distribution func- 
tion of the type similar to Eq. ( |2.9D . In the classical limit 
the density-matrix equation goes over to the FPE and, 
accordingly, the quantum LLB simplifies to Eq. ( 2.16| ) 
with the microscopically determined bath-coupling pa- 
rameter A. For ferromagnets with the dominant homo- 
geneous exchange interaction the quantu m LL B equation 
simplifies to the same macroscopic form ( 4.13 ). The main 
result of the present paper — the "slow" LLB equation 
( |3.17 ) — can be obtained in the quantum case, too, by a 
perturbative solution of the density-matrix equation for 
a slowly changing magnetic field, which is similar to the 
derivation in Sec. [II. The final result can be, however, 
obtained by replacing Fi ^ Fl int = in the longitudi- 
nal relaxation term. The analytical expression for Fi^int 
in the quantum case without single-site anisotropy was 
already given in Ref. ^. Very recently it was generalized 
for the anisotropic case to describe the-jthermoactivation 
escape rate of quantum spin systems.t3 

The most serious problem by the derivation of the LLB 
equation for ferromagnets, that has not been solved yet, 
is taking into account the spin-spin interactions. This is 
a rather difficult task, since one should consider the FPE 



(4.5) for the whole system, which describes all possible 
static and dynamic spin correlations. Even at an equi- 
librium, where th e so lution of the FPE (4.5) is known 



and given by Eq. ( |4.q ) , one faces the problem of a phase 
transition in a many-body system. Calculation of spin- 
spin contributions into the longitudinal and transverse 
kinetic coefficients Li a nd L2 in the LLB equation for 
ferromagnets, Eq. (4.13), goes beyond the scope of this 
paper and is planned for the future. 
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APPENDIX: DERIVATION OF THE 
FOKKER-PLANCK EQUATION 



Here the derivation of the Fokker-Planck equation (2.2 ) 
is presented, which isjiiore direct and simple than the 
original one by Browned and which uses more advanced 
stochastic methods applied, in particulafmin the dynami- 
cal renormalization-group (RG) theoryE'Era The RG con- 
siderations start, however, with "soft-spin" models with 
the formal Langevin sources (i.e., the inhomogeneous 
terms in the stochastic differential equations for the spin 
density), which cannot be interpreted as random mag- 
netic fields acting on spins. For our purposes, we will 
derive the FPE for magnetic systems with the methods 
of Refs. ^ but starting from the more realistic stochas- 
tic Landau-Lifshitz equation (2T). At first we introduce 
the probability distribution of the random Gaussian noise 

c, 



■ exp 



1 

2a 



drC^r) 



(Al) 



where = JD^T is the noise partition function, JDC 
denotes functional integration over realizations of Ci''') 
and a = 2XT/{jfio). With the help of Eq. (^ the 
average of any noise functional ^[C] can be written as 



With the use of the obvious identity 



SapS{T - t) 



one can calculate variations of of Eq. ( |Al[ ): 



(5Ca(i) a 



^Ca(i)^C/3(i') 

etc. Since for all n one has 



DC 



0, 



(A2) 
(A3) 
(A4) 

(A5) 



the functio nal integration of Eq. (A4) leads to {Ca{t)) = 
and Eq. ( |2.2| ). Further, one can show that all averages 
of an odd number of ^ components are zero and those 
of an even number n > 2 of C's decay pairwise and can 



be expressed through the pair average Eq. (2.2), i.e., the 
statistics of the random field C{t) is Gaussian. 

The distribution function of spins / is determined as 

/(N,i) = (7r(t,[C]))c, 7T{t,[C])=S{N-s{t)). (A6) 

The time derivative of / can be calculated using 



dN 



(A7) 



and the equation of motion (2.1), which yields 



a/ 

dt 



d 



7[N X H]/ - 7A[N X [N X H]]/ 



-7[Nx (C(i)7r(i, [C]))< 



(A8) 



Then the average {C,{t)n{t, \C\))c can be transformed with 
the use of the first of Eqs. ( A_4) and integration by parts. 



(C(t)^(t, [C]))< 



-aJoCnit, [C]) 



a{ - 



\ Kit) 



6C{t) 
I dTi 6s pit, [C]) 



(A9) 



where with a = x,y,z are the orts of the Descarte 
coordinate system and summation over components a, (3 
is implied. The variational derivative 6sf3/6Ca can be 
calculated, if one writes down the for ma l solution of the 
stochastic Landau-Lifshitz equation (O), 



0it)=j dt'ep^a.s^it')[H,,it') + Cc.it')] + 



(AlO) 



where ep^a is the antisymmetric unit tensor. One can 
see that 



6Ca.it') 



jepjaS^it'), t' <t 

t' >t. 



(All) 



For t ~t' the above calculation does not yield a definite 
value of 6s[3/6Ca, but with the help of the usual argu- 
ments based on the regularization of 6 functiongJ the lat- 
ter can be found to be ■^jeap-fSj{t). Now Eq. (A9) can 
be finally written in the form 



(C(t)^(t, [C]))c = Y 



N X 



or 

ON 



(A12) 



Adopting it in Eq. (|A8|), one comes to the Fokker-Planck 
equation (2^). 
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